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A Ward-Takahashi identity preserving Bethe-Salpeter kernel can always be calculated explicitly 
from a dressed-quark-gluon vertex whose diagrammatic content is enumerable. We illustrate that 
fact using a vertex obtained via the complete resummation of dressed-gluon ladders. While this 
vertex is planar, the vertex-consistent kernel is nonplanar and that is true for any dressed vertex. 
In an exemplifying model the rainbow-ladder truncation of the gap and Bethe-Salpeter equations 
yields many results; e.g., 7T- and p-meson masses, that are changed little by including higher-order 
corrections. Repulsion generated by nonplanar diagrams in the vertex-consistent Bethe-Salpeter 
kernel for quark-quark scattering is sufficient to guarantee that diquark bound states do not exist. 



PACS numbers: 12.38.Aw, 11.30.Rd, 12.38.Lg, 24.85. +p 

I. INTRODUCTION 

Dynamical chiral symmetry breaking (DCSB) and con- 
finement are keystones in an understanding of strong in- 
teraction observables and their explanation via a nonper- 
turbative treatment of QCD. The gap equation Q] 
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is an insightful tool that has long been used to explore 
the connection between these phenomena and the long- 
range behaviour of the interaction in QCD ||. In this 
equation: D^ LV (k) is the renormalised dressed-gluon prop- 
agator; T®(q;p) is the renormalised dressed-quark-gluon 
vertex; mt,m is the A-dependent current-quark bare mass 
that appears in the Lagrangian; and :— J A d 4 g/(27r) 4 
represents a translationally-invariant regularisation of 
the integral, with A the regularisation mass-scale. The 
quark-gluon-vertex and quark wave function renormali- 
sation constants, Zi(£ 2 ,A 2 ) and Z2(£ 2 ,A 2 ) respectively, 
depend on the renormalisation point and the regularisa- 
tion mass-scale. 

The solution of Eq. ([!]) is the dressed-quark propaga- 
tor, which takes the form 
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and is obtained by solving the gap equation subject to the 
renormalisation condition that at some large, spacelike C 2 



S(p) 1 L, 2=C 2 = H ■ P + m(Q , 



(3) 



where m(C) is the renormalised current-quark mass at the 
scale (: Z±m(£) = Zi mb m , with Z4 the renormalisation 



constant for the scalar part of the quark self energy. At 
one loop in perturbation theory 



m(C) 
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where: rh is the renormalisation-group-invariant current- 
quark mass; j m = 12/(33 — 2 Nf) is the leading-order 
mass anomalous dimension, for Nf active flavours; and 
Aqcd is the A/-flavour QCD mass-scale. 

Since QCD is an asymptotically free theory the chiral 
limit is unambiguously defined by rh — ||, which can 
be implemented in Eq. (|l|) by applying [Q 



lim Z 2 (C 2 ,A 2 )m(A) = 0. 



(5) 



The formation of a gap, described by Eq. (|l]), is identified 
with the appearance of a solution for the dressed-quark 
propagator in which m(C) ~ 0(1/C 2 ) 7^ 0; i.e., a solu- 
tion in which the mass function is power-law suppressed. 
This is DCSB. It is impossible at any finite order in per- 
turbation theory and entails the appearance of a nonzero 
value for the vacuum quark condensate: 

,A 

-(qq)° c = lim Z 4 (C 2 , A 2 ) A c tr D / S°(q,0, (6) 

A^oo J q 

where tr£> identifies a trace over Dirac indices alone and 
the superscript "0" indicates the quantity was calculated 
in the chiral limit. 

It is apparent that the kernel of the gap equation is 
formed from a product of the dressed-gluon propagator 
and dressed-quark-gluon vertex. The kernel may be cal- 
culated in perturbation theory but that is inadequate for 
the study of intrinsically nonperturbative phenomena. 
Consequently, to make model-independent statements 
about DCSB one must employ an alternative systematic 
and chiral symmetry preserving truncation scheme. 
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One such scheme was introduced in Ref. ||. Its 
leading-order term is the rainbow-ladder truncation of 
the DSEs and the general procedure provides a means to 
identify, a priori, those channels in which that truncation 
is likely to be accurate. This scheme underlies the suc- 
cessful application of a renormalisation-group-improved 
rainbow-ladder model to flavour-nonsinglet pseudoscalar 
mesons |Q and vector mesons ||, fj], |[ g], and indicates 
why the leading-order truncation is inadequate for scalar 
mesons and flavour-singlet pseudoscalars [|l0). The sys- 
tematic nature of the scheme has also made possible a 
proof of Goldstone's theorem in QCD ||. 

In quantitative applications, however, the leading- 
order term alone has been used almost exclusively: 
Refs. ^ are exceptions but they consider just the 
next-to-leading-order term. Hence one goal of our study 
is a nonperturbative verification of the leading-order 
truncation's accuracy. 

One element of the gap equation's kernel is the dressed- 
gluon propagator, which in Landau gauge can be written 



J fJ,V 
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It has been the focus of DSE studies 
simulations 
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Jl2j and lattice 
and contemporary analyses suggest 



that d(k 2 X 2 )/k 2 is finite, and of 0(1 GeV" 2 ), at k 2 = 0. 
However, this behaviour is difficult to reconcile with the 
existence and mag nitude of DCSB in the strong interac- 
tion spectrum ]la| : it is a model- independent result that 
a description of observable phenomena requires a kernel 
in the gap equation with significant integrated strength 
on the domain k 2 < 1 GeV 2 |l6| . The required magnifica- 
tion may arise via an enhancement in the dressed-quark- 
gluon vertex but, hitherto, no calculation of the vertex 
exhibits such behaviour JT^ ] . Hence another aim of our 
study is to contribute to the store of nonperturbative in- 
formation about this vertex. 

In Sec. |n| we briefly recapitulate on the truncation 
scheme of Ref. ||. Then, using a simple confining 
model [Ul], we demonstrate that an infinite subclass 
of contributions to the dressed-quark-gluon vertex: the 
dressed-gluon ladders, can be resummed via an algebraic 
recursion relation, which provides a closed form result for 
the vertex expressed solely in terms of the dressed-quark 
propagator. This facilitates a simultaneous solution of 
the coupled gap and vertex equations obtained via the 
infinite rcsummation, as we describe in Sec. III. While 



the algebraic simplicity of these results is peculiar to our 
rudimentary model, we anticipate that the qualitative 
behaviour of the solutions is not. 

In Sec. [V we describe the general procedure that 
enables a calculation of the Bethe-Salpeter kernel for 
flavour-nonsinglet mesons that is consistent with the 
fully-resummed dressed-gluon-ladder vertex. The kernel 
is itself a resummation of infinitely many diagrams; and it 
is not planar, an outcome necessary to ensure the preser- 
vation of Ward-Takahashi identities. This kernel is the 
heart of the inhomogeneous vertex equations and asso- 



ciated bound state equations whose solutions relate to 
strong interaction observables. In our simple model it, 
like the vertex, can also be obtained via an algebraic re- 
cursion relation, complete and in a practical closed form. 
In Sec. ^ we also study the bound state equations in 
a number of meson channels, and derive and solve the 
analogous equation for diquark channels. 
Section [v| is a summary. 



II. A DRESSED-QUARK-GLUON VERTEX 

The truncation scheme introduced in Ref. j^] may be 
described as a dressed-loop expansion of the dressed- 
quark-gluon vertices that appear in the half ampu- 
tated dressed-quark-antiquark (or -quark-quark) scat- 
tering matrix: S 2 K, a renormalisation-group invari- 
ant, where K is the dressed-quark-antiquark (or -quark- 
quark) scattering kernel. All n-point functions involved 
thereafter in connecting two particular quark-gluon ver- 
tices are fully dressed. The effect of this truncation in 
the gap equation, Eq. ([!]), is realised through the fol- 
lowing representation of the dressed-quark-gluon vertex, 
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Here T 3g is the dressed-three-gluon vertex and it is ap- 
parent that the lowest order contribution to each term 
written explicitly is 0(g 2 ). The ellipsis represents terms 
whose leading contribution is 0(<? 4 ); i.e., the crossed-box 
and two-rung dressed-gluon ladder diagrams, and also 
terms of higher leading-order. 

The expansion of S 2 K just described, with its implica- 
tions for other n-point functions; e.g., the dressed-quark- 
photon vertex, yields an ordered truncation of the DSEs 
that, term- by-term, guarantees the preservation of vec- 
tor and axial- vector Ward-Takahashi identities, a feature 
exploited in Ref. || to prove Goldstone's theorem. Fur- 
thermore, it is readily seen that inserting Eq. (||) into 
Eq. (0) provides the rule by which the renormalisation- 
group-improved rainbow-ladder truncation |^, ^, ^, [i| 
can be systematically improved. It thereby facilitates an 
explicit enumeration of corrections to the impulse current 
that is widely used in calculations of clcctrowcak hadron 
form factors H . 



A. Resumming Dressed-Gluon Ladders 

We can't say anything about a complete resummation 
of the terms in Eq. (p|). However, we are able to con- 
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FIG. 1: Integral equation for a planar dressed-quark-gluon 
vertex obtained by neglecting contributions associated with 
explicit gluon self-interactions. Solid circles indicate fully 
dressed propagators. The vertices are not dressed. 



tribute to aspects of the more modest problem obtained 
by retaining only the sum of dressed-gluon ladders; i.e., 
aspects of the vertex depicted in Fig. [IJ This infinite 
subclass of diagrams is l/7V c -suppressed. 



since, in general, one can write 

r n(p) = "i(p 2 )7m + a 2 (p 2 )j -pp^ 

-a 3 (p 2 )ip p + a 4 (p 2 )«7 M 7 -p, (10) 

but it does not restrict our ability to address the ques- 
tions we raised in the introduction because those am- 
plitudes which survive are the most significant in the 
dressed-quark-photon vertex and it is an enhancement 
in the vicinity of p 2 = that may be important for a re- 
alisation of DCSB using an infrared-finite dressed-gluon 
propagator Jl5[ . 

The summation depicted in Fig. [l] is expressed via 

r M (p+,p-) = ZfSiU 

+ 1 J V pa {p- Z) lp S{l + )T^l + ,l_)S{l_) la 

(11) 



1. A Model 

To simplify our analysis and make the key elements 
transparent we employ the confining model introduced 
in Ref. |l8|, which is defined by the following choice for 
the dressed-gluon line in Fig. m. 



2V(fc) :=<T2V(fc) = <W 



(9) 



Plainly, Q, measured in GeV, sets the model's mass-scale 
and henceforth we set Q = 1 so that all mass-dimensioned 
quantities are measured in units of Q. In the following, 
since the model is ultraviolet-finite, we usually remove 
the regularisation mass-scale to infinity and set the renor- 
malisation constants equal to one. 

The model defined by Eq. (|9|) is a precursor to an ef- 
ficacious class of models that employ a renormalisation- 
group-improved effective interaction and whose contem- 
porary application is reviewed in Refs. |ll| |2(J. It has 
many positive features in common with that class and, 
furthermore, its particular momentum-dependence works 
to advantage in reducing integral equations to algebraic 
equations that preserve the character of the original equa- 
tion. Naturally, there is a drawback: the simple mo- 
mentum dependence also leads to some mo del- dependent 
artefacts, but they are easily identified and hence are not 
cause for concern. 



2. Planar Vertex 

The general form of the dressed-quark gluon vertex 
involves twelve distinct scalar form factors but using 
Eq. (0) that part of this vertex which contributes to the 
gap equation has no dependence on the total momentum 
of the quark-antiquark pair; i.e., only T p (p) := T p {p,p) 
contributes. This considerably simplifies the analysis 



with v± — vi-P, where v = £,p, etc., is any four-vector, 
but using Eq. (|]) this simplifies to 

r P (p) = 7p + 1 i P s(p) r M (p) s( P ) lp , (12) 

where the additional factor of 3/4 (1/8 cf. 1/6) owes 
itself to the combined operation of the (5-function and 
the longitudinal projection operator. Inserting Eq. ( |To| ) 
into Eq. (|l^) one finds «4 = and hence the solution of 
Eq. (O) simplifies: 



r /x(p) = "i(p 2 )7m + a 2(p 2 )-f -PPix - a3(p 2 )iPn- (13) 
We re-express this vertex as 
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(14) 



i=0 



Ki(p 2 )7m + ^{V 2 )! -PPp ~ OL\{p 2 )ipy\ 



(15) 



where the superscript enumerates the order of the iterate: 
r^ =0 is the bare vertex, 



a\ = 1 , a° 2 = - a° 3 ; 



(16) 



is the result of inserting this into the r.h.s. of Eq (|12j) 
to obtain the one-rung dressed-gluon correction; T^ =2 is 
the result of inserting and is therefore the two-rung 
dressed-gluon correction; etc. 

Now a simple but important observation is that each 
iterate is related to its precursor via the following recur- 
sion relation: 

r7 +1 b + ,P-) :=( a r' +1 ( P+lP _) = 



2V(p - u) S(u+) l a T p (u+, u _) S(u_) l b lu , 

(17) 
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FIG. 2: Recursion relation for the iterates in the fully- 
resummed dressed-gluon-ladder vertex, Eq. (|l7|): filled circles 
denote a fully-dressed propagator or vertex! Using Eq. 
p = k, and this relation is expressed by Eq. (|li|). 



which is depicted in Fig. g. Using Eq. (^) this simplifies: 

r p +1 (p) = 1 -^s( P )r p (p)s(p) 1 , (is) 

and substituting Eq. ( |l5| ) into Eq. ( |l8| ) yields (s = p 2 ) 



ai +1 ( S ) 

a l+1 (s):=| a l +\ s ) | = C(s; A, B) a l (s) 



(19) 



with this vertex. Using Eq. @ that gap equation is 

s(p) _ 1= (int-P + m + KSWTM (23) 

I i7-p + m + r jLl (p)S'(p)7 A1 



and substituting Eq. (|l3|) this gives 
1 



sA 2 + B 2 
1 

sA 2 + B 2 



[A(2a 1 -sa 2 )-Ba s ] , 



(24) 



[B (4ai + sa2) — sA as] 



(25) 



A(s) = 1 + 
B(s) = rnd 



Obviously, Eqs. (B4J), (E5|), completed using Eqs. (|2£ 
form a closed, algebraic system. It can easily be solved 
numerically and that procedure yields simultaneously 
the complete gluon-ladder-dressed vertex and the prop- 
agator for a quark fully dressed via gluons coupling 
through this nonpcrturbative vertex. 

We note here that in the chiral limit a realisation of 
chiral symmetry in the Wigner-Weyl mode is always pos- 
sible. This realisation is expressed via the B = solu- 
tion of the gap equation, which from Eqs. (p2|), (25) is 
evidently always admissible for m = 0. 



0(s;A,B) = 
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\ AAB AsAB 2(B 2 - sA 2 ) 



Now it is clear that Eqs. @, (p|), @ entail 
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(20) 



(21) 



where the last step is valid whenever an iterative solution 
of Eq. ( |Tl] ) exists, and defines a solution otherwise, so 
that, with A = sA 2 (s) + B 2 (s), 



Oil 
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4A 
1 + 4A ' 

-8 A 2 1 + 2 A 

l + 2(B 2 -,sA 2 )-8A 2 1 + 4 A' 

-8AB 
1 + 2{B 2 -sA 2 )-8A 2 ' 



(22) 



We have thus arrived at a closed form for the 
gluon-ladder-dressed quark-gluon vertex of Fig. [l| i.e., 
Eqs. (|l3|), (|22|). Its momentum-dependence is determined 
by that of the dressed-quark propagator, which is ob- 
tained by solving the gap equation, itself constructed 



III. SOLUTIONS OF THE GAP AND VERTEX 
EQUATIONS 

A. Algebraic Results 

Before reporting results obtained via a numerical solu- 
tion we consider a special case that signals the magnitude 
of the effects produced by the complete resummation in 
Fig. [|; i.e., we focus on the solutions at s = 0. In this 
instance Eqs. © give, with A = A(0), B = B(0), 



AB 2 



1 



A B r 



8 A 2 



2B 2 



1 



ai(s = 0) 
a 2 {s = 0) 
a 3 (s = 0) 

Substituting these expressions, Eq. ( |25| ) becomes 

16 Bo 



(AB 2 + l) 2 2B 2 -1 
8A a B a 1 
AB 2 + 1 2B 2 -1' 



(26) 



B 



AB 2 



1 



(27) 



and in the chiral limit this yields 

B = -VTE^ 1.94, 



(28) 



which makes plain that the model specified by Eq. (||) 
supports a realisation of chiral symmetry in the Nambu- 
Goldstone mode; i.e., DCSB. The value in Eq. (p8[) can 
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be compared with that obtained using the bare vertex; 
i.e., the leading-order term in the truncation of Ref. ||: 



B, 



(0) 




2, to see that the completely resummed dressed- 



gluon-ladder vertex alters Bq by only 
Similarly, Eq. (|24|) becomes 



-3%. 



A = 1 



8A 



4B 2 



1 



1 



1 



2B 2 -1 



which using Eq. (E8|) gives 



A 



26 
15 



1.73. 



(29) 



(30) 



This, too, may be compared with the leading-order re- 
2. Again the resummation does not materi- 



sult: A 



(0) 




ally affect the value: here the change is 
Inserting Eqs. ©, (M) into Eqs. (f26| 



■13%. 
one finds 



«i(0) 
a 2 (0) 
« 3 (0) 



15 

16 
221 
1800 
1 

%/l5 



0.94 cf. 1.0. 
0.12 cf. 0.0. 
0.26 cf. 0.0. 



(31) 



where the last, comparative column lists the values for 
the leading-order (bare) vertex. It is evident that the 
solution of Eq. ([ll]) obtained using an infrared amplified 
effective interaction, Eq. (^), which supports DCSB (and 
also confinement (2^]), does not exhibit an enhancement 
in a neighbourhood of p 2 = 0. This provides an internally 
consistent picture: a dressed-quark-antiquark scattering 
kernel, whose embedding in the gap equation already pos- 
sesses sufficient integrable strength, does not additionally 
magnify itself. 
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FIG. 3: Numerical solution of Eq. (^2|): solid line; cf. the 
solution obtained with the bare vertex: dotted line; and the 
one-loop corrected vertex: dashed line. 



which has two solutions: ka 
(physical) solution of Eq. (|3 
Therefore sA 2 (s) 



sA 2 



3/4, 5/4. The required 
satisfies A(s) s zl^° i+ 
from above and hence the 



physical branch of the solution is described by 



K A 



< 2. 



(35) 



Evidently, in the chirally symmetric case too, the ladder- 
dressed vertex reduces the magnitude of the solution. 

These features are apparent in the complete numerical 
solution of Eq. (|j2"|), which is depicted in Fi g, p . Further- 
more, the gluon-ladder-dressed vertex, Eq. (|13|), obtained 
with this solution does not exhibit an enhancement. 



2. Nambu-Goldstone Mode 



B. Numerical Results 



Wigner-Weyl Mode 



The B = solution of Eq. (E5() is always admitted 
when m = 0. In that case Eq. (BJ) becomes 



A(s) 



1 



8A(s) 



l + AsA 2 {s) 



1 



1 



1 -AsA 2 {s) 



(32) 



of which there is no closed-form solution. (It is a quintic 
equation for A(s).) 

However, using the bare vertex one finds a chiral limit 
solution: A^ \s) — \/2/s for s ps 0. We therefore sup- 
pose that in the neighbourhood of s = Eq. ( |32"| ) admits 
a solution of the form 



A. 



-(*)=m 1/2 . 



V s / 

Substituting this into Eq. ( |32"| ) yields 
16 1~2k a 



1 



1 + 4 ka 1 - 4 ka 



(33) 



(34) 



In Sec. HI A we described features of the Nambu- 



Goldstone mode solution of the gap equation, whose 
p 2 ~ properties are smoothly related to those of the 
m 7^ solution, as is also the case in QCD (^2). A 
complete solution is only available numerically and our 
calculated results for the dressed-quark propagator and 
gluon-ladder-dressed vertex are depicted in Figs. || - [| 

It is apparent from the figures that the complete re- 
summation of dressed-gluon ladders yields a result for 
the dressed-quark propagator that is little different from 
that obtained with the one-loop-corrected vertex; and 
there is no material difference from the result obtained 
using the zeroth-order vertex. A single, exemplary quan- 
tification of this observation is provided by a comparison 
between the values of M(s — 0) = B(0)/A(0) calculated 
using vertices dressed at different orders: 



M(0) 



N = N = 1 A = 2 A = oo 



1 



1.105 1.115 1.117 



(36) 



Similar observations apply to the vertex itself. Of 
course, there is a qualitative difference between the 
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FIG. 4: Upper panel - A(s), M(s) obtained with m = 0: 
solid line; Lower panel - A(s), M(s) obtained with m — 0.023: 
solid line. All dimensioned quantities are expressed in units of 
Q in Eq. (S). A fit to meson observables requires Q ~ 0.5 GeV 
and hence this value of the current-quark mass corresponds 
to ~ 10 MeV. In both panels, for comparison, we also plot the 
results obtained with the zeroth-order vertex: dotted line; and 
the one-loop vertex: dashed line. 
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FIG. 5: The dressed-quark propagator can be written: 
S(p) — —i r y-pav(s) + as{s). Here we plot 4 ov(s), 2crs(s) 
obtained with m — 0.023 cf. the results obtained with the 
zeroth-order vertex: dotted line; and the one-loop vertex: 
dashed line. Clearly, the analytic properties of the dressed- 
quark propagator are qualitatively unaffected by our dressing 
of the vertex. 
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zeroth-order vertex and the one- loop-corrected result: 
a 2.3 7^ in the latter case. However, once that effect 
is seeded, the higher-loop corrections do little. 



FIG. 6: aj, i = 1,2,3, calculated with m = 0.023. These 
functions calculated at one-loop (dotted line) and two-loop 
(dashed line) are also plotted for comparison. The results 
obtained with m = are little different. 



3. Vertex Ansatz 

In the absence of a nonperturbatively dressed quark- 
gluon vertex a number of phenomenological DSE studies 
have employed an Ansatz, which is based on a nonpertur- 
bative analysis of QED and, in particular, is constrained 
by the vector Ward-Takahashi identity and the require- 
ment of multiplicative renormalisability |p3| . That vertex 
has been used J24| with the model interaction of Eq. (|^) 
in which its form is expressed via 

/ \ \ , s n dA(s) , . n dB(s) 

oi(a) = 4 s , oa(fl) = 2 , a 3 (s) = 2 — . 

as as 

(37) 

The functions defined by these expressions, calculated us- 
ing the self-consistent solutions determined in Ref. 
are depicted in Fig. They bear little resemblance to 



the functions obtained systematically via the resumma- 
tion of dressed-gluon ladders, an outcome that could be 
anticipated based on the difference between the dressed- 
quark propagator calculated herein and that in Ref. [M . 
This finding does not invalidate the Ansatz nor its use in 
modelling QCD but merely shows that the Ansatz cannot 
be a sum of dressed-gluon ladders alone. 



IV. BETHE-SALPETER EQUATION 

The renormalised homogeneous Bethe-Salpeter equa- 
tion (BSE) for the quark-antiquark channel denoted by 
M can compactly be expressed as 

,A 

[r M (k;P)] E F= [K(k,q;P)]^[ X M(q;P)]GH (38) 
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FIG. 7: Qj, i = 1,2,3, calculated from Eqs. ( |37| ) using the 
dressed-quark propagator solutions obtained in Ref. 



where: TM{k;P) is the meson's Bethe-Salpeter ampli- 
tude, k is the relative momentum of the quark-antiquark 
pair and P is their total momentum; E, . . . .H represent 
colour, flavour and spinor indices; and 



XM(k;P) = S(k+)T M (k;P) S(fc_) . 



(39) 



In Eq. (|3q) , which is depicted in Fig. |§|, K is the 
fully-amputated dressed-quark-antiquark scattering ker- 
nel. The choice 



[K(k, q; P)]%f = V^k - q) [l a ^] EG [l a lu 



1HF i 



(40) 



yields the dressed-gluon ladder-truncation of the BSE, 
which provides the foundation for many contemporary, 
field-theory-based phenomenological studies of meson 
properties, see, e.g., Refs. p5fl. 



A. Vertex-consistent Kernel 

The preservation of Ward-Takahashi identities in those 
channels related to strong interaction observables re- 
quires a conspiracy between the dressed-quark-gluon ver- 
tex and the Bethe-Salpeter kernel ]^[ ^6). We now de- 
scribe a systematic procedure for building that kernel. 

As described, e.g., in Ref. gfj], the DSE for the dressed- 
quark propagator, S, is expressed via 



ST[S] 
SS 



= 0, 



(41) 



where T[S] is a Cornwall- Jackiw-Tomboulis- like effective 
action. The Bethe-Salpeter kernel is then obtained via 
an additional functional derivative: 



k gh 

^EF 



5T, 



EF 



SS, 



GH 



(42) 



Herein the self energy is given by the gap equation, 
Eq. (|2^), and the recursive nature of the dressed-gluon- 
ladder vertex entails that the n-th order contribution to 





FIG. 8: Homogeneous BSE, Eq. (|3J). Filled circles: dressed 
propagators or vertices; K is the dressed-quark-antiquark 
scattering kernel. A systematic truncation of S 2 K is the key 
to preserving Ward-Takahashi identities [H, G6|. 



the kernel is obtained from the rt-loop contribution to the 
self energy: 



D F (j,-?)!°T,%)i a r;(?,p). (43) 



Since T M (p, q) itself depends on S then Eq. ( f42| ) yields a 
sum of two terms: 

[K n (k,q;P)]<i« 
= 2V(*-9)P°7,Jeg [l a rS(q-,k-)] HF 



s 



EL 



SS GH {q±) 



[i°r:(L,L)] 



LF 



(44) 



Here, in addition to the usual effect of differentiation, the 
functional derivative adds P to the argument of every 
quark line through which it is commuted when applying 
the product rule. NB. T)^ also depends on S because of 
quark vacuum polarisation diagrams. However, as noted 
in Ref. || , the additional term arising from the derivative 
of T>^ v does not contribute to the BSE kernel for flavour 
non-diagonal systems, which are our focus herein, and 
hence is neglected for simplicity. It must be included 
though to obtain a kernel adequate for an analysis of 
problems such as the 77-77' mass splitting, for example. 
Now, introducing [xAffe P)]gh, the BSE becomes 



T M (k-P)= f V^(k-q) 

J a 



X*°7/. 



XM (q;P)l a T v (g- , fc_ ) + S(q + ) k a Mv (q,k;P) , 

(45) 



where we have used Eq. (jl4|) and defined 

A a Mv ( q> k;P) = J2^l(<l>k;P), (46) 



n=0 



with 



[AS(^;P)] 



LF 



q SS(q±) GH 



[l a Tl{l_,k_)] LF [ X M(q;P)} 



\GH ■ 



(47) 
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Equation ( |45| ) is depicted in Fig. The first term 
is instantly available once one has an explicit form for 
T". To develop an understanding of the second term, 
which is identified by the shaded box in the figure, we 



employ the recursive expression for the dressed-quark- 
gluon vertex, Eq. JTij), with p — k, P = I — k, to obtain 
an inhomogeneous recursion relation for A^™ : 



J 



A%l(i,k;P) = / V p<T (£~q)l b lp XM(q;P)l a T:- 1 (q-,q-+k-e)S(q-+k-£)l b la 



+ / V pa (k - q) l h lp S(q+ +£-k) l a K-\q + + I - k, q + ) X mW, P) 1*1* 

l-A 

V pa (£ - q')l h lp S(q' + ) A a AT\q', q' + k-l-P) S{q'_ + k - t) l b Ja . 



r 



(48) 



Equation (Eq) is illustrated in Fig. [l^ and, combined 
with Figs. §, this discloses the content of the vertex- 
consistent Bethe-Salpeter kernel: namely, it consists of a 
countable infinity of contributions, an infinite subclass of 
which are crossed-ladder diagrams and hence nonplanar. 
It is clear that every n > 0, vertex-consistent kernel must 
contain nonplanar diagrams. Charge conjugation can be 
used to expose a diagrammatic symmetry in the Bethe- 
Salpeter kernel, which is the procedure we used, e.g., to 
obtain Eq. (p3|). (The steps and outcome described here 
formalise the procedure illustrated in Fig. f of Ref. f|.) 

At this point if T™ and the propagator functions: A, B, 
are known, then A\j v , and hence the channel-projected 
Bethe-Salpeter kernel, can be calculated explicitly. That 
needs to be done separately for each channel because, 
e.g., Kiu depends on xm{t,P)- 

To proceed we observe that the Bethe-Salpeter ampli- 
tude for a 7r-meson is 



TUk;P)=I c ri l5 [ifl(k 2 ,k-P;P 2 ) 

+ 1 -Pfl{k\k-P-P 2 )+ 1 -kk-Pfl{k\k-P;P 2 ) 
+ a plJ k p P u f^(k 2 ,k-P;P 2 )] , (49) 

where I c is the identity in colour space (mesons are colour 
singlets) and {t j ; j — 1,2,3} are the Pauli matrices. This 
illustrates the general structure of meson Bethe-Salpeter 
amplitudes, which hereafter we express via 

N M 

T M (k; P) = IcT G l (k; P) P M {k 2 , k ■ P; P 2 ) , (50) 



where Q l (k;P) are the independent Dirac matrices re- 
quired to span the space containing the meson under 
consideration. (Subsequently, to simplify our analysis, 
we focus on Nf — 2 and assume isospin symmetry. There 
is no impediment, in principle, to a generalisation.) 

The result of substituting Eq. ( |50|) into the r.h.s. of 
Eq. (Eq) can be expressed in the compact form 



f = Hf=(H 1 +H 2 )f; 



(51) 



where / is a column vector composed of the scalar func- 
tions: f i = f % M . Here Hi is the contribution from the 
first term on the r.h.s. in Eq. ( Eq 1 ) and it is a Nm x Nm 
matrix wherein the elements of row-j are obtained via 
that Dirac trace projection which yields f M on the l.h.s.; 
i.e., for Vi such that 



fi = tTD [VjTm] 



then, using l a l a = -C 2 {R) = -| for SU(N C = 3), 



(52) 



V pv (k-£) lp ,S{£+) 



x g k (t;P)S(t-)T v {t-,k-) 



f k M (£ 2 A-P;P z ) 



(53) 



H.2 represents the contribution from the second term, 
to which we now turn. In mesonic channels the colour 
structure of A 1 ^,™ is simple: 



lv Mu — 1 lv Mu 



(54) 







FIG. 9: BSE expressed in Eq. (]4J), which is valid whenever FIG. 10: Recursion relation for A^™, Eq. (|4J). We label the 

diagrams on the r.h.s., from left to right, as Qi, Q2, C. 



Y p can be obtained via a recursion relation 
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because xm oc I c and l b l a l b = l a /(2N c ). The other part 
of the direct product is a matrix in Dirac space that can 
be decomposed as follows: 

N A 

ASr„ = J2(%&k;P)g*(£,k;P) (55) 

A=l 

N A N m 
A=l j=l 

where the sum over j implicitly expresses an integral 
over the relative momentum appearing in \m {gi A = 
1, . . . , 2Va} are the independent Dirac matrices required 
to completely describe A^, whose form and number 
are determined by the structure of XMj > Nm> an d 
{/3™(P 2 ); A = 1, . . . , A^a} are the associated scalar coef- 
ficient functions. (We subsequently suppress momentum 
arguments and integrations for notational ease.) 

Using Eq. (|5^) , the recursion relation of Eq. pq) trans- 
lates into a relation for {/3™; i = A, . . . , A^a}. To obtain 
that relation one first isolates these functions via trace 
projections. That can be achieved by using any com- 
plete set of projection operators: {P\. v \ i = 1, • • ■ , N\}, 
in which case one has 



where O resolves the dressed-quark-gluon-vertex recur- 
sion, as we saw with our simple dressed-gluon model in 
connection with Eq. ([l9]), and the first term vanishes be- 
cause (3° = by definition, Eq. @. 

Finally the complete BSE involves the sum expressed 
in Eq. (Eq), which is determined by: 

oo oc n — 1 

(3 = ^2 /3 n = t Mi ] J MGO n -i- x oP 

n— 1 71=1 j=0 



i=0 



1 



1 ML 



MG 



1-0 



a 



(62) 



and this, via Eq (p5|), completely determines the second 
term in Eq. (Ea) so that we can complete Eq. d5lh with 



x ^S(l + )g*(e,k;P)0*(l,k;P) 



f h Ai(P 2 ), (63) 



(3 n x = [M] X \> tr cd 



A 



Mu 



(57) 



where trc^ identifies a trace over colour and Dirac in- 
dices. (NB. The optimal choice of projection opera- 
tors would yield [.M]aa' = 5\\<, as we assumed, e.g., 
in Eq. (|52"|).) We subsequently adopt a compact matrix 
representation of Eq. (|57]) : 



(3 n = MT n 



(58) 



where /3 n is a column vector with A^a entries. 

Replacing A^™ on the r.h.s. in Eq. ([57]) by the r.h.s. of 
Eq. ([l8|) and using the distributive property of the trace 
operation, one obtains 



T 



a 



Lf3 n - 



(59) 



which we have now demonstrated is calculable in a closed 
form. (Recall that the sum over k implicitly expresses 
an integral over the relative momentum in the Bethe- 
Salpeter amplitude. Note, too, that there is no sum over 
iterates in this equation: it is (3 from Eq. (|6^) which 
appears.) 



B. Solutions of the vertex-consistent meson 
Bethe-Salpeter Equation 

To elucidate the content of the BSE just derived we re- 
turn to the algebraic model generated by Eq. (^) . In this 
case solutions of the BSE are required to have relative 
momentum k = so that Eq. (|5p| ) simplifies to 



where G describes the contributions to the trace from 
the first two terms in Fig. |To|, G\ t 2, which are determined 
by the dressed-quark-gluon vertex and thus proportional 
to a™ -1 , and L represents the contribution from the last 
term, C. Using Eq. (p9|), Eq. (p8|) becomes 



(60) 



(3 n = M [G a' 1 - 1 + L (3 n - 1 } 



(We reiterate that M, G and L are all functionals of xm 
and hence are different for each meson channel.) 
It is evident that Eq. (160) entails 



n-l 

f3 n = [ML] n (3° + E[MIf MG a™^ 1 

3=0 



= [ ML Y MGO n - j - l a° , 

3=0 



(61) 



r 



M 



N<N M 

(P)= 5 l (P)h(P 2 ) (64) 



and Eq. (|5l| ) is truly algebraic; i.e., there are no implicit 
integrations. Furthermore, the kernel H = H(P 2 ); i.e., 
it is a matrix valued function of P 2 alone, and therefore 
the mass, M%, of any bound state solution is determined 
by the condition 



det [H(P 2 ) - 1] 



P 2 +Mjj=0 



0. 



(65) 



which is the requirement for any matrix equation: Hx = 
x, to have a nontrivial solution: it is the characteristic 
equation. NB. If no solution of Eq. ( p5| ) exists then the 
model doesn't produce a bound state in the channel un- 
der consideration. 
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Tv-meson 



In our algebraic model, because k = 0, Eq. 
plifies to 

r w (P)=75 \ih(P 2 )+l-Pf2(P 2 



sim- 



(66) 



where P is the direction-vector associated with P; P 2 = 
1, and for the projection operators of Eq. ([32]) we choose 



Pi 



■75 , P2 



■ 7 ' P 75 • 



(67) 



The vertex A° ; ™ transforms as an axial-vector and its 
form is therefore spanned by twelve independent Dirac 
amplitudes. However, since the relative momentum is 
required to vanish, that simplifies and we have: 

A?" = Z°75 [ff(P 2 )~/ u + (%(P 2 ) 7 -PP v 

+ (P 2 ) P v + PI (P 2 ) 7 „ 7 • P] , (68) 

and an obvious choice for the projection operators of 
Eq. @ is 



pi 



- 2 1 P,7-^75 



7^75, PL ;I / 



P,75 



(69) 



1- Plul5- 



It can now be shown that G^ = in Eq. ( |59| ) : is the 
sum of two terms, the first and second in Eq. (p8|), and 
using Eq. (^) and recalling that this model forces the rel- 
ative momentum to vanish in the bound state amplitude, 
these two terms are equal in magnitude but opposite in 
sign because the projection operators are axial-vector in 
character. (At leading order this result expresses an ex- 
act cancellation between one of the one-loop vertex cor- 
rections and the crossed-box term in Fig. 1 of Ref. ||.) 
Gtt = is a consequence of the fact that under charge 
conjugation A 7r/J (£, k;P) transforms according to: 



A^(-fc,-^;P) t --A 7rAI (£,A:;P), 



(70) 



where (•) denotes matrix transpose, coupled with the 
result that Eq. (§) enforces k = = I in the BSE. It 
is not a general feature of the vertex-consistent Bethe- 
Salpeter kernel. 

It is plain from Eqs. (pa), (p2n that G T = entails 



A a = 



(71) 



and hence the complete vertex-consistent pion BSE is 
simply / = Hif; i.e., 

TAP) = - 7m S(Q) TAP) S(-Q) TA-Q, -Q) , (72) 

with Q = P/2, and the dressed-quark propagator and 
dressed-quark gluon vertex calculated in Sec. 1IIB 2. 



The characteristic polynomial obtained from Eq. (|72j) 
is plotted in Fig. [ll]: the zero gives the pion's mass, 



1j& 








Bare 




1 loop 




2 loop 




Recursive 



FIG. 11: Characteristic polynomial obtained in the chi- 
ral limit from the vertex-consistent BSEs for the it- and p- 
mesons: the pion is plainly massless. The function obtained 
using m = 0.023, which corresponds to w 10 MeV, is almost 
indistinguishable on the scale of this figure. 



Eq. ([35|) , which is listed in Table |. Figure [IT] provides a 
forthright demonstration that the pion is massless in the 
chiral limit. This is a model- independent consequence 
of the consistency between the Bethe-Salpeter kernel we 
have constructed and the kernel in the gap equation. 
The figure illustrates that the vertex-consistent Bethe- 
Salpeter kernel converges just as rapidly as the dressed- 
vertex itself (cf. Fig. |TJ). The insensitivity, evident in the 
table, of the pion's mass to the order of the truncation is 
also explained by the fact that our construction preserves 
the axial-vector Ward-Takahashi identity. 



2. p-meson 

The complete form of the Bethe-Salpeter amplitude for 
a vector meson in our algebraic model is 



^(P)= 7 -£ A (P)/T(P 2 



a, v eAP)P u f^P 2 ). (73) 



This expression, which only has two independent func- 
tions, is much simpler than that allowed by a more re- 
alistic interaction, wherein there are eight terms. Nev- 
ertheless, Eq. ( |7q ) retains the amplitudes that, in more 



TABLE I: Calculated n and p meson masses, in GeV, quoted 
with Q = 0.48 GeV, in which case m = 0.023 5 = 11 MeV. (In 
the notation of Ref. Js), this value of g corresponds to 77 = 
0.96 GeV.) n is the number of dressed-gluon rungs retained in 
the planar vertex, see Fig. [lj and hence the order of the vertex- 
consistent Bethe-Salpeter kernel: the rapid convergence of the 
kernel is apparent in the tabulated results. 





M%=° 


M£ =1 


M£= 2 


MJT 00 


7T, 


m 


= 














7T, 


in 


= 0.011 


0.152 


0.152 


0.152 


0.152 


P, 


m 


= 


0.678 


0.745 


0.754 


0.754 


(>■ 


m 


= 0.011 


0.695 


0.762 


0.770 


0.770 
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sophisticated studies, are found to be dominant |pj. In 
Eq. @, {e A (P); A = -1,0, +1} is the polarisation four- 
vector: 



P-e A (P) = 0, VA; e A (P) • e A ' (P) = S 



rAA' 



Here the projection operators of Eq. (52) are 



P A = ± 7 • z\P) > = ~ e A (P) P 
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(74) 



(75) 



To solve the vertex-consistent BSE we need to calculate 
A p ; "(e A , P), which is defined in Eq. (|48|) and depends on 
the Bethe-Salpeter amplitude in Eq. fl73). We have from 
Eq. (||) that 



A^( e A ,P) = l a A n pu (e\P), 



(76) 



and in our algebraic model the Dirac structure is com- 
pletely expressed through 

+ %(P) * 7 • e x P v + ^ 4 (P) a aP e x a (P) Pp 7 „ 

+ /? p " 5 (P) io af} e x a (P) P P P v + /3 p " 6 (P) 7 • e A 7 „ , (77) 

which has this simple form because A pl ,(e A ,P) cannot 
depend on the relative momenta. In this case an obvious 
choice for the projection operators of Eq. (|57|) is 



'A -7 



P5 

' A.;i/ 



4 



P 2 



- 7 • Pe, 

4 ' I - 



P.7-^ 



■7, £a(P)P/3, 

-Ip.^e^P)^, P? ;I/ 



(78) 



The substitution of Eqs. @, @ into Eq. @ yields 



M = - 

p 2 



3 




1 3 -1 
-1 -1 1 



1-1 

3-1-1 

1 

1 



-1 3 
-1 1 



(79) 



which simply expresses the fact that, e.g. 
1 



P1(P) 



tr cd 



1 A p ;w 



Pa 



(80) 



The next step is a determination of the matrix G p 
in Eq. (|59|), which gives the contribution to the kernel's 
recursion relation from the first two terms in Fig. [Tc| . 
That is achieved by substituting Eq. (|76|) into Eq. (48), 



and using Eq. ^ and its consequences this yields 
2A(Q 2 )M P G p = 

2{f?(P)B(Q 2 ) + tfQA(Q 2 )} 

fi(P) Q MQ 2 ) - fZ{ p ) b(Q 2 ) 

-f?(P)QA(Q 2 ) + f?(P)B(Q 2 
-tf(P)QA(Q 2 ) + f%(P)B(Q 2 





(81) 



with Q = VQ 2 = v / P 2 /2, and A(s) = sA 2 (s) + B 2 (s) as 
in Eq. (p^). It is immediately apparent that the Dirac 
components associated with /3 p5 , /3 p6 are annihilated by 
this part of the interaction. Furthermore, the complete 
contribution to /3™ is 



M P G P 



-\Q 2 ) 
-\Q 2 ) 
-\Q 2 ) 



(82) 



and so it is evident that the subleading Dirac compo- 
nents of the dressed-quark-gluon vertex do not contribute 
to the A°-related part of the vertex-consistent Bethe- 
Salpeter kernel: they are eliminated by the last two 
columns in M p G p . (These two simplifications are a fea- 
ture of our model.) 

The final element we require is the matrix L p in 



Eq. (59); i.e., the contribution to the kernel from the 
last term in Fig. |Io|. That is obtained by substituting 
Eq. (j7q) into the last term of Eq. ( |48"| ) , which gives 

4 A 2 M p L p = 



- 2A 














2A 





-A 





-2A 














2[Q 2 A 2 - P 2 ] 


2P 2 


-2QAB 


-2QAB 











A 













































(83) 

where the argument of each function is Q 2 ; e.g., B = 
B(Q 2 ). In our representation M p L p does not exhibit an 
explicit dependence on r p (P). That dependence is ac- 
quired through the recursion relation since (3 p cx M p G p , 
as is apparent from Eq. (|60[). 

Note that this part of the kernel also annihilates the 
Dirac components associated with /3 p5 , /3 p6 in A p;1 ,(A, P). 
Hence the active form of the vertex is 

A^( e \ P) = /^(P)e A + /5» 2 (P) etij-P 

+ (3%{P) i 1 -e x P v + /?« 4 (P) a afi e x (P) P p 7 „, (84) 

where the /3™, i — 1, . . . ,4 are obtained from Eq. (|6£ 
using the elements calculated above. 

Putting all this together we have the complete vertex- 
consistent BSE for the p-meson: 

T A (P) = -^S(Q)T X (P)S(~Q)T^~Q,-Q) 

- 7 ^(Q)A pp (e A ;P), (85) 
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with the dressed-quark propa gator a nd dressed-quark 
gluon vertex calculated in Sec. IIIB 2 , and A pp (0,0; P) 
obtained from Eqs. @, (||), (|62|) using the matrices 
displayed above. 

The characteristic polynomial obtained from Eq. (|85| ) 
is plotted in Fig. Its zero gives the p-meson mass and 
that is listed in Table |. The tabulated values demon- 
strate that the rainbow-ladder truncation underestimates 
the result of the complete calculation by < 10% (i.e., the 
calculation performed using the fully-resummed dressed- 
gluon-ladder vertex) and simply including the consis- 
tent one-loop corrections to the quark-gluon-vertex and 
Bethe-Salpeter kernel reduces that discrepancy to < 1%- 
Furthermore, it is clear that at every order of trunca- 
tion the bulk of m p is obtained in the chiral limit, which 
emphasises that the ir-p mass-splitting is driven by the 
DCSB mechanism. 



3. Scalar and Axial-vector Mesons 

In a simple constituent-quark picture, the ground state 
scalar and axial vector mesons are angular-momentum 
L = 1 eigen states. This qualitative feature is ex- 
pressed in their Poincare covariant Bethe-Salpeter am- 
plitudes through the presence of materially important 
relative-momentum-dependent Dirac components; e.g., 
Refs. [ p7| . However, the model defined by Eq. (||) forces 
meson Bethe-Salpeter amplitudes to be independent of 
the constituent's relative momentum and, owing primar- 
ily to that, the rainbow-ladder truncation of the model 
generates neither scalar nor axial-vector meson bound 
states Jig . We have found, unsurprisingly, that improv- 
ing the description of the dressed-quark-gluon vertex and 
Bethe-Salpeter kernel is insufficient to overcome this de- 
fect of the model. However, we anticipate that these 
improvements will materially alter the results of BSE 
studies of scalar mesons that employ more realistic in- 
teractions; e.g., Refs. 1 28 . 



C. Diquark Bethe-Salpeter Equation 

The ladder-rainbow truncation, which is obtained by 
keeping only the first term in Eq. (^|), generates colour- 
antitriplet quark-quark (diquark) bound states [^9| . Such 
states are not observed in the hadron spectrum, and it 
was demonstrated in Ref. || that they are not present 
when one employs the one-loop-dressed vertex and con- 
sistent Bethe-Salpeter kernel. Herein we can verify that 
this feature persists with the complete planar vertex and 
consistent kernel. To continue, however, we must slightly 
modify the procedure described above because of the 
colour-antitriplet nature of the diquark correlation. (NB. 
Colour-sextet states are not bound in any truncation be- 
cause even single gluon exchange is repulsive in this chan- 
nel. It is for this reason, too, that colour-octet mesons 
do not appear. Note also that the absence of colour- 



antitriplet diquark bound states does not preclude the 
possibility that correlations in this channel may play an 
important role in nucleon structure fl30[ since some at- 
traction does exist, e.g., Ref. pi]] .) 
The analogue of Eq. (|3^) is 

r A - 

[r D (k:P)] EF = [K(k,q;P)]^[ XD (q;P)] G H (86) 

Jq 

where: k is the relative momentum of the quark-quark 
pair and P is their total momentum, as before; and 



X D(q;P) = S(q + )T D (q:P) 



(87) 



with Tr)(q; P) the putative diquark's Bethe-Salpeter am- 
plitude. In this case K is the fully-amputated dressed- 
quark-quark-scattering kernel, for which 



K(k,q;P)]%f = 2V(* ~ <z) Tl^EG [(l a l»? 



HF 



(88) 



yields the dressed-gluon ladder-truncation of the BSE. 

Inspection of the general structure of K reveals that, 
following our ordering of diagrams, it can be obtained 
directly from the kernel in the meson BSE via the re- 
placement 



S(fc)i°7M- [7/* 1° 



(89) 



in each antiquark segment of K, which can be traced 
unambiguously from the external antiquark line of the 
meson's Bethe-Salpeter amplitude. 

The appearance here of a matrix transpose makes the 
definition of a modified Bethe-Salpeter amplitude p9[] 
useful: 



T c D (k-P) :=T D {k-P) & 



(90) 



where C — 7274 is the charge conjugation matrix, from 
which it follows that 



x g(fc;P) - S{k+)T%{k;P)S(k-), 



(91) 



using C7*C t = — 7 M . rg(fe;P) satisfies a BSE whose 
Dirac structure is identical to that of the meson BSE. 
However, its colour structure is different, with a factor 
of (— replacing l a at every gluon vertex on what was 
the conjugate-quark leg. (NB. It is from this modification 
that the absence of diquark bound states must arise using 
the dressed-ladder vertex, if it arises at all.) For example, 
if one has a term in the meson BSE of the form 

l a ^ S l h lv S l c lp S T M Sl a lp S l c lp S l b 7 u (92) 

then the related term in the equation for r£(fc; P) is 

i°7 M s i b lv s i c lp sv c D s (-/°)*7/i s {-rfi P s {-l b f lv . 

(93) 

There are three colour-antitriplet diquarks and their 
colour structure is described by the matrices 

{A^ ; k = 1, 2, 3 : = A 7 , A 2 A = A 5 , A 3 A = A 2 } , (94) 
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and, as for mesons, their Bethe-Salpeter amplitudes can the rainbow-ladder meson BSE except that 
be written as a direct product, colour ® Dirac: 



rg(fc;P) = A*r£(fc;P). 



(95) 



l a I c l a = -^I c 



l a A* (-l a f 



- 3 ^a- (96) 



For colour-singlet mesons the colour factor is simply the 

identity matrix I c . (Recall that we focus on Nf = 2 and In both cases the colour matrix now factorises and can 

assume isospin symmetry. Hence the diquark's flavour therefore be cancelled. Hence the rainbow-ladder BSEs 

structure, which is described by the Pauli matrix r 2 in satisfied by the colour- independent parts of Tm and 

this case, cancels in the BSE.) r£(fc;P) are identical but for a 50% reduction of the 

As we observed, the rainbow-ladder diquark BSE is coupling in the diquark equation. This expresses the 

obtained by using Eq. ( {88| ) in Eq. (|S6|). Right-multiplying fact that ladder-like dressed-gluon exchange between two 

the equation thus obtained by C'we find immediately quarks is attractive and explains the existence of diquark 

that the equation satisfied by Tp(k;P) is the same as bound states in this truncation p9| . 

I 



Following the above discussion it is apparent that the diquark counterpart of Eq. ( J48| ) is 

A a £(e, k- p) = tv pa {i - q ) i b lp x %(q; p) rrV, <?- + k-e) (i a Y s( Q - + k-t) (i b f la 

J a 



V pa {k - q) l b lp S(q+ +£-k) l a T n v -\q + + t - k, q+) X %(q; P) (HS 



V p *(t - q')l b l P S(q' + ) L (q\ q' + k - £; P) S(q'_ + k - I) (l b )H 



(97) 



The factorisation of the colour structure observed in the ladder-like diquark BSE persists at higher orders and can be 
used to obtain a recursion relation for the diquark kernel analogous to that depicted in Fig. 10. Indeed, a consideration 
of Eq. (|9^) reveals that in general one can write 

A% n v (£, k; P) = AUi, k; P) l a X A + A$ v (l, k; P) X k A (l a f , (98) 

where in this direct product A™ v and AJ, are Dirac matrices. Defining: 



7'L:n \k ja A a ' n TOA n ah 1 \fi 



V prT (£ - q) lp X C qq {q\ P) rr \q-, Q- + k - t) S(q- +k-l) la 



V pa {k - q) 7p S(q+ +t-k) K~ 1 (q + + t - k, q+) P) 7 , 



V pa {i-q') lp S{q') ~±A'?- 1 (q',q' + k-£;P)-?-A%- 1 (q',q' + k-e;P) S(q'_ + k - I) la \ A*, (99) 



5 a n — 1 / / / 



T** X k := AZ (O* = - ~ l~A n lv + 2 A»J A* 



v P „{t - q ) i P Xaaiq; p) rr 1 (?- ,g- + k-£) + k - 1) 7ct 



+ °- / V pa (k~q) lp S(q + +£-k)T:- i (q + +e-k,q + )x^ q (q;P)l, 



+ I V pa {£-q') lp S(q' + ) 'lA^iq^q' + k-tP)- L -A n 2 -\q\q l + k-i-P) S{q'_ + k - i) la \ A*, (100) 



1 \n—lf„> J 



then 



NB. The first lines in each of Eqs. (99), (J100|) prove that 



An 1 \r-fh\n | n rrR'n] \n 1 \r\r-j-'L;n , s-r-R;n~\ 

\v — — 2 l J v Z 1 v J ' iV 2u — ~ 2 iv J • 

(101) 
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there is no mixing between colour-antitriplet and colour- 
sextet diquarks, whose colour structure is described by 
the six symmetric Gell-Mann matrices. 

We can now write the vertex-consistent BSE for the 



colour-antitriplet diquark channels (cf. Eq. (45)) 



X7^ [x^;P)K(£-,k_) - S(£ + )T v L ' n (£,k;P)] . 

(102) 

It is straightforward to verify that this equation repro- 
duces the diagrams considered explicitly in Ref. |l| . 

Having factorised the colour str uctu re, the summation 
appearing in the first term of Eq. (102) yields the dressed 
vertex we have already calculated. One proceeds with the 
second term by analogy with Eq. ( |55| ) and observes that 
the matrix-valued functions, A12, can be decomposed: 



A?„= E r% x {t,k;P)h x v {t,k;P)- 



(103) 



where {h* ; A = 1, . . . , Na d } is the smallest set of Dirac 
matrices capable of expressing A^J™ completely, whose 
form and number depend on the channel under consider- 
ation. Projection operators, Vp. v , are easily constructed 
so that, with 



77 := column(?/ii, 
we have 



,TllN An , V21, • ■ ■ ,r)2N An ), (104) 



[ V n ] A = [M'] AA dr D 7>*A?„ 



(105) 



where A, A' = 1, . . . , 2N Ad , and I = 1 for A' < N A and 
I = 2 for A' > N\ D . Now the procedure of Eqs. (|58| ) - 
( |60| ) can be repeated to arrive at 



r) n = M' [G' a"- 1 + L' rj"- 1 ] 



(106) 



where G' , L' are natural analogues of the matrices in- 
troduced in Eq. (|59|), and thereafter one can continue to 
obtain an obvious extension of Eq. ( J62j). T his completely 
determines the second term in Eq. ( |102| ) and thus we 
have arrived at the vertex-consistent BSE for the colour- 
antitriplet diquark channel. 



D. Solutions of the Diquark Equation 

1. Scalar Diquark 

In the algebraic model specified by Eq. (||) the general 
Dirac structure of a J p — + quark-quark correlation is 



r° + (P) = 75 ff{P 2 ) + i-Pff{P 2 ) 



(107) 



\ \ \\ 

\ N \ 








Bare 

1 loop 

2 loop 

Recursive 



FIG. 12: The characteristic polynomial obtained using 
m = 0.023, which corresponds to ~ lOMeV, calculated for 
the scalar and axial- vector co lour -antitriplet diquark channels 
using the diquark BSE, Eq. (Fo|). 



which is the same as that of the pion, for reasons which 
are obvious given the discussion following Eq. (po|). The 
characteristic equation for this channel is obtain ed fol- 
lowing the method made explicit in Sec. IV Bl and it 
is depicted in Fig. [12]. The existence of a bound state 
for n — 0; i.e., in the rainbow-ladder truncation, is ap- 
parent. However, so too is the effect of the higher-order 
terms, which was identified in Ref. Q: at each higher- 
order nonplanar diagrams in the kernel provide signifi- 
cant repulsion, which overwhelms any attraction at that 
and preceding orders and thereby ensures diquark con- 
finement; i.e., the absence of coloured quark-quark bound 
states in the spectrum. This feature is retained by the 
completely resummed kernel, using which, instead of a 
zero, the characteristic polynomial exhibits a pole: the 
repulsion is consummated. (This feature is not tied to 
the interaction in Eq. (S); e.g., Ref. [p2|.) 



2. Axial Vector Diquark 

The general form of a J p = 1 + quark-quark correlation 
is represented by 

rJJ(P) = 7 ■ t x (P) ff(P 2 ) + 4(P) Pu ft(P 2 ) ■ 

(108) 

Calculating the characteristic polynomial for this channel 
is also a straightforward application of methods already 
introduced and it is plotted in Fig. |lj. The features in 
this channel are qualitatively identical to those of the 
scalar diquark. 



V. SUMMARY 

Using a planar quark-gluon vertex obtained through 
the resummation of dressed-gluon ladders we have ex- 
plicitly demonstrated that from a dressed-quark-gluon 
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vertex, obtained via an enumerable series of terms, it 
is always possible to construct a vertex-consistent Bethe- 
Salpeter kernel that ensures the preservation of Ward- 
Takahashi identities in the physical channels related to 
strong interaction observables. While we employed a 
rudimentary model to make the construction transpar- 
ent, the procedure is general. However, the algebraic 
simplicity of the analysis is peculiar to our model. For 
example, using a more realistic interaction the gap and 
vertex equations would yield a system of twelve coupled 
integral equations. Nevertheless, we anticipate that the 
qualitative features highlighted herein arc robust. 

The simple interaction we employed characterises a 
class of models in which the kernel of the gap equa- 
tion has sufficient integrated strength to support dy- 
namical chiral symmetry breaking (DCSB). The com- 
plete ladder summation of this interaction, calculated 
self-consistently with the solution of the gap equation, 
produces a dressed-vertex that is little changed cf. the 
bare vertex. In particular, it does not exhibit an enhance- 
ment in the vicinity of k 2 = 0, where k is the momen- 
tum carried by the model dressed-gluon. In addition, the 
dressed-quark propagator obtained in this self-consistent 
solution is qualitatively indistinguishable from that ob- 
tained using the rainbow truncation. 

The vertex-consistent Bethe-Salpeter kernel is neces- 
sarily nonplanar, even when the vertex itself is planar, 
and in our simple model it is easily calculable in a closed 
form: for a more realistic interaction it can be obtained 
as the solution of a determined integral equation. The 
fact that our construction ensures the kernel's consis- 
tency with the vertex and hence preservation of Ward- 
Takahashi identities is manifest in the Goldstone boson 
nature of pion, which is preserved order-by-order and in 
the infinite resummation. 

Our explicit calculations focused primarily on flavour- 



nonsinglet pseudoscalar mesons and vector mesons. We 
found that a consistent, nonperturbative dressing of the 
vertex and kernel changes the masses of these mesons by 
< 10% cf. the values obtained using the rainbow-ladder 
truncation. That is not the case in the pseudoscalar 
channel if the kernel is dressed inconsistently. Further- 
more, 90% of the w-p mass-splitting is already generated 
in the rainbow-ladder truncation, which emphasises that 
this splitting is primarily driven by DCSB. The rainbow- 
ladder truncation is a poor approximation for flavour- 
singlet pseudoscalar mesons and scalar mesons. 

We also considered quark-quark scattering and found 
that, with anything but a ladder-like vertex-consistent 
Bethe-Salpeter kernel, diquark bound states do not exist 
in the spectrum. 
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